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Abstract 

Recent XMM-Newton observations of clusters of galaxies have indicated the soft X-ray spectra to be 
inconsistent with the simple isobaric cooling flow model. There is almost no feature of the cold gas expected 
from the model. This shows that we have not yet understood the physics of the hot gas in clusters of galaxies 
well. A quantitative evaluation of the behavior of gas cooling is important not only for understanding the 
clusters, themselves, but also for studying cosmology and galaxy formation. To clarify the problem of this 
reported discrepancy, we have studied scaling relations for clusters of galaxies based on the self-similarity 
assumption. We also propose an observational strategy to solve this problem. 

Key words: galaxies:clusters:general — galaxiesxlustersxooling flows — cosmologydarge-scale struc- 
ture of universe 



1. Introduction 

Previous soft X-ray imaging studies with medium 
spectral resolution have shown evidence of cool- 
ing flows in many clusters of galaxies. In many 
cases, a simple isobaric cooling flow model (|^00 
[cite] cite.jfet92 Johnstone et al. (1992)) was used, 
and showed a reasonable fit to the observed re- 
sults. However, recent XMM-Newton observations 
with high spectral resolution (^00 [citeJcite.tamOlTamura 
et al. (2001); ^0 [citcJcite.kaaOlKaastra ct al. 
(2001); 000 [cite]cite.pet01Peterson et al. (2001); §00 
[cite] cite. mbfbOlMatsushita et al. (2001)) showed that 
the soft X-ray spectra is inconsistent with the simple iso- 
baric cooling flow model. In A1795, Tamura et al. (2001) 
showed that the XMM-Newton RGS spectrum obtained 
from the cluster center can be described by an isother- 
mal model with kT ~ 4 keV, and that the upper limit 
of the emission measure of the cool component {kT < 1 
keV) is only a few %. Kaastra et al. (2001) showed a 
significant lack of cool gas below 1.5 keV in the Sersic 
159-03 cluster of galaxies, where a cooling flow of 230 Mq 
yr~^ has been indicated from the ROSAT PSPC obser- 
vation (^0 [cite]cite.all97Allen, Fabian (1997)). Peterson 
et al. (2001) showed that the RGS energy spectrum at 
the cluster center of A1835 requires a cut-off in the emis- 
sion measure distribution at 2.7 keV to obtain a good fit 
with the isobaric cooling flow model. In the Virgo cluster, 
Matsushita et al. (2001) reported that the emission mea- 
sure of the cool component below 1 keV is much lower than 
the expected value from the cooling flow model. In sum- 
mary, although the central excess of the surface brightness 
profile has suggested cooling flow and the existence of cold 
gas, RGS high-resolution spectra at the center of cooling 
flow clusters did not exhibit the expected emission lines 



from the cold gas, and the temperature distribution had 
to be cut-off (cut-off temperature, Tcut) at 1-3 kcV (here- 
after, we call this feature the "lack of cold gas" problem). 
This indicates that we have not yet understood the physics 
of the hot gas in clusters well. 

A proper understanding of the process of gas cooling is 
important in many respects. First, it is necessary to dis- 
cuss the cluster gas evolution and its status (for example, 
whether the hot gas is in hydrostatic equilibrium or not). 
Second, it is a key to understand the dark matter prop- 
erties and evolution, since we cannot observe dark matter 
directly, and its information is obtained only through the 
observed hot gas properties. Third, clusters are impor- 
tant to study cosmology (for example, to estimate the 
cosmological parameters), and have been widely used in 
research. In order to use clusters as a cosmological probe, 
we need to understand the cluster properties well. Fourth, 
a study of the cooling gas in clusters gives important infor- 
mation for understanding galaxy formation and evolution, 
because we can regard galaxies as small mass clusters as 
a first approximation when we discuss galaxy formation 
and evolution. Then, cooling is more important in galax- 
ies than in clusters, since the cooling time of the hot gas 
in galaxy-size halos is shorter than the cosmic time (this 
feature is different from the clusters' case) and stars are 
formed from the cooled gas (we have obtained various in- 
formation about galaxies through stars). 

Thus, it is important to understand the cluster 
hot gas properties. In particular, we need to solve 
this "lack of cold gas" problem. There are sev- 
eral possibilities to solve this problem: heating, ther- 
mal conduction, inhomogencous mctallicity, absorption 
by cold gas, and non-standard cooling function, and 
so on (0)0 [citc]citc.fab01Fabian et al. (2001); 000 
[citeJcite.petOlPeterson et al. (2001)). These mechanisms 
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have both advantages and disadvantages, and it is difficult 
to judge which is the main mechanism at present. 

It seems to be difficult to solve the "lack of cold gas" 
problem by studying individual clusters, as in previous 
literature. Thus, in this Letter, we discuss this problem 
while considering correlations between cluster properties. 

2. Scaling Relation 

In this section, we consider correlations between cluster 
properties, particularly between Lcooi (as defined below) 
and the gas temperature (or cluster mass). First, we re- 
view the self-similarity model briefly and then modify it 
in order to explain the observed luminosity-temperature 
relation. We then consider the cooling effect on the cluster 
properties. 

At first, we consider the whole properties of clus- 
ters. If shock heating caused by gravitational collapse 
is the dominant mechanism, it is thus naturally ex- 
pected that clusters follow a simple scaling model (^0 
[cite]cite.kai86Kaiser (1986)). We assume that clusters 
of galaxies are spherically symmetric, and that the hot 
gas is isothermal and in a steady state with hydrostatic 
equilibrium. We also assume that the gas number den- 
sity profile n(r) is described by the conventional (3 model, 
n(r) = no/[l + {r /r^fff^/'^ , with /3 = 2/3, and that the 
core radius is proportional to the cluster radius RcV- 
Tc oc i?ci- Assuming that the average mass density in clus- 
ters and the gas fraction are identical in all clusters, we 
obtain the following relations: 



TT-o oc (= constant), Tc oc M^Y^, 



2/3 



where Md and Tvir are the cluster mass and virial tem- 
perature, respectively. Hereafter, we assume that the gas 
temperature, Tgas, equals the virial temperature (Tgas = 
Tvir) when gas cooling is neglected. The above relations 

1/2 

imply that Tc oc Rd oc Tgas • Although whether the obser- 
vational data satisfies this relation is still in controversy; 
for example, Vikhlinin et al. (1999) found that the radius 

1/2 

of a fixed mean gas overdensity of 1000 is scaled to Tgas • 
If brcmsstrahlung emission is the dominant emis- 
sion mechanism, the bolometric luminosity, Tboi, is 
proportional to the gas temperature squared: Tboi oc 



no Tc" J gas OC Tgj^g. However, the observed luminosity- 
temperature relation is steeper than this relation; the lu- 
minosity is roughly proportional to the gas temperature 
cubed: iboi oc T'^^^ (e.g., |o|DO [cite]cite.dsj93David et al. 
(1993)). There arc many arguments to explain this dis- 
crepancy. We can divide them into roughly two groups. 
One is that the gas fraction varies with the cluster mass; 
the other is that the gas profile deviates from the self- 
similarity [equation (1)]. For example, the core radius is 
not proportional to i?ci, or the cluster radius is not propor- 
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tional to Tgas ■ We briefly discuss the second possibility 
in section 4. Here, we take the first possibility. We aban- 
don the assumption of the gas fraction being constant, 
and allow it to vary as a function of the cluster mass, so 
that the observed luminosity-temperature relation is re- 



produced (e.g., ^0 [cite]cite.dsj93David et al. (1993)). 
This is expected if non-gravitational processes affect hot 
gas properties. The pre-heating model is an example (see 
e.g., ^0 [cite]citc.bcm01Bialek et al. (2001)), and the 
variable galaxy formation efficiency is another example 
(see e.g., |^00 [cite] cite. ptceOOPearce et al. (2000); ^0 
[cite] cite. mtkOlMuanwong et al. (2001)). If the galaxy 
formation efficiency is more efficient in poor clusters than 
in rich systems, the gas fraction increases with the clus- 
ter mass. In order to reproduce the observed luminosity- 
temperature relation, the gas fraction must be propor- 



tional to Tga4' {Mll^y. 



no oc , 



0CM„V^ TvirOcAf^^ 



(2) 



cl ' ^vir ^ J"cl 

We use these relations in the following discussion. 

Next, we consider the cooling effect. We assumed the 
hot gas to be isothermal in the above discussion. However, 
at the center, the cooling time may be shorter than the 
cluster age and the cooling mechanism becomes impor- 
tant. We define the cooling radius, rcooi, as the radius 
where the cooling time, TcooI, equals the cluster age, Tago: 
■'■cooi(''cooi) ~ Tagc- Here, the cooling time is defined as 
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where e = An^Tgas'' is the bolometric emissivity and A is a 
numerical constant. Then, the cooling radius is described 

as 

1/2 
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(4) 



It is to be noted that the value in the bracket does not 

1 /2 

depend on the cluster mass, since no oc Tg4 . Then, the 
gas number density at TcooI is given by 



'^(fcool) 
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[1 + (rcooi/rc)2 
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(5) 



Hereafter, Tgas denotes the gas temperature outside the 
cooling radius. At r < TcooI, the gas profile may devi- 
ate from the /3-model. As a fiducial model, we take the 
self-similar gas density and temperature profile within the 
cooling radius. 



T(r) =Tgas/T 



^ ^cool 

(for r < Tcooi), 



n{r) = n{r cool) fn 



\ ^cool 



(6) 



where /t and /„ are functions independent of cluster 
mass. We then obtain a relation between the luminos- 
ity within the cooling radius, Lcooi, and the outer-region 
gas temperature, Tgas (or the cluster mass Md), as 

/•I-cool 

icooi= / Airr^ e {r) dr (xn{r cool frl^oiTg^ 

(7) 

That is, the ratio of Tcooi to Tboi is constant. In this case, 
the cut-off temperature, Tcut, must be scaled to Tgas; then. 



oc T3,, cx Ml 



No. ] 

it is expected that the ratio Tcut/?gas is the same in all 
clusters. 

At present, there is no compilation of the correlation 
between icooi defined above (or Lcooi/iboi) and Tgas- 
Therefore, in order to give a rough estimate, we use the 
data compiled by Allen (2000, table 8) here. It is to be 
noted that icooi, Alien is estimated from the cooling flow 
model. Lcooi,Aiion differs from Lcooi, defined above in gen- 
eral, but we expect that there is a positive correlation 
between Lcooi,Aiicn and icooi- His results showed that 
there is no clear trend between -/jcooi,Aiicn/^boi and Tgas, 
although the scatter is large. This may indicate that clus- 
ters are roughly described by the scaling relation [equation 
(6)]. 

3. Heating Models 

The above discussion provides useful information for 
approaching the "lack of the cold gas" problem. In this 
section, as an example, we consider the possibility that 
heating is the main mechanism to solve the problem. The 
heating model has one possible advantage in that it nat- 
urally explains that there is no, or little, mass deposition 
at the cluster center observationally. Of course, there are 
problems; for example, the required energy is huge. If we 
assume a heating source of AkT kcV in the cluster core 
with a radius of rioox 100 kpc and a gas number density 
of na X 10~^ cm~^ in rg Gyr, the required average energy 
input rate is 

U _ SnAkTV 

Since detailed modeling is beyond the scope of this Letter, 
we do not discuss this problem any more. 

Because cooling flow features are seen in many clusters, 
it is reasonable to assume that they arc in a steady state. 
Within the cooling radius, the cooling time is shorter than 
the cluster age, by definition. In order that clusters are 
in a steady state, the cooling has to be balanced with 
heating. If clusters are self-similar within the cooling ra- 
dius [equation (6)], then the heating rate has to be pro- 
portional to the gas temperature cubed (the cluster mass 
squared) from the result obtained in the previous section. 
Hereafter, we assume that the heating rate, F, is described 
by r cx Affj using a parameter, p. 

When the exponent of p equals to 2, clusters can become 
self-similar within the cooling radius [equation (6)] . In this 
case, the "lack of cold gas" problem have to be observed 
in all cooling flow clusters and the ratio Tcut/Jgas is the 
same in every cooling flow cluster. 

In the case of p < 2, the heating affects the hot gas more 
significantly in poor clusters than in rich clusters. Thus, 
in poor clusters, the ratio of the cut-off temperature to 
the gas temperature outside the cooling radius becomes 
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larger, that is, the temperature decrement to the center 
becomes weaker. In rich clusters, the simple isobaric cool- 
ing flow model may be a relatively good approximation. 

On the contrary, in the case of p > 2, the heating signif- 
icantly affects the hot gas in rich clusters. Then, in poor 
clusters the "lack of cold gas" problem disappears and the 
simple isobaric cooling flow model is applicable. 

At present, several candidates for the heating source 
are proposed. The above discussion can tell what kind 
of properties are required in the heating sources. The 
most popular candidate is supernova. If we consider only 
normal star formation, the supernova rate is roughly pro- 
portional to the total cluster optical luminosity. Lci, or 
the cluster mass: the heating rate due to supernova is pro- 
portional to the cluster mass, that is p = 1. In this case, 
the "lack of cold gas" problem should be observed only in 
poor clusters. 

Active Galactic Nuclei activity is another candidate. 
The total optical luminosity of a cD galaxy Lcd is propor- 
tional to Lji^s ^0 [cite]cite.bahOOBahcall (2000)). If the 
activity of cD galaxies (e.g., jet kinetic energy) is propor- 
tional to the optical luminosity, Lcd , then the heating rate 
follows p^l. There is another possibility. The cD galax- 
ies are more extended than other giant elliptical galaxies; 
they show very extended low surface brightness envelopes. 
The luminosity of the cD galaxy's envelope, Lenv, is pro- 
portional to (|qo [cite]cite.bahOOBahcall (2000)). If 
heating is an accompanying result of the formation of the 
cD galaxy's envelope, then the exponent becomes p ~ 2. 
Then, clusters can become self-similar within the cooling 
radius [equation (6)]. 

Next, we consider galaxy-galaxy collisions. If the ki- 
netic energy of galaxies is released in galaxy-galaxy colli- 
sions, r AgaingaiO'ci(mgai(T^i/2), whcrc Agai is the num- 
ber of galaxies in a cluster and Cci is the velocity disper- 

1 /3 

sion of the galaxies. Since A^gai cx Md and CTci csc , 
the heating rate, F, is ^ M^^; that is, p = 2. In this case, 
clusters can become self-similar within the cooling radius 
[equation (6)]. 

Finally, we consider heating by the magnetohydrody- 
namic effects proposed by Makishima (see e.g., Makishima 
et al. 2001). The galaxy motion causes significant mag- 
netohydrodynamic turbulence and frequent magnetic re- 
connection. Since the kinetic energy of the galaxies is the 
origin of the heating, the heating rate may be proportional 
to iVgaiCTci- Then, the heating rate, F, is (x A^gaicr^i/rheating, 
where Thoating is the heating time scale. Since we do not 
know the detail of the heating, we consider the following 
two simple cases of the heating time scale. One is that the 
heating time scale is proportional to the crossing time, i.e., 
Rci/o'ci- The other case is that the heating time scale is 
constant. In both cases, the heating rate, F, is propor- 

5 /3 

tional to . Then, the "lack of cold gas" problem may 
be observed in most clusters; the problem rather tends to 
become weak in poor clusters. 

We only consider cooling flow clusters in the above dis- 
cussion. However, some candidates of the heating process 
(e.g., supernovae and galaxy-galaxy collisions) also take 
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place in the non-cooling flow clusters. They may affect 
the hot gas properties (such as the temperature structure) 
in the non-cooling flow clusters as well as in the cooling 
flow clusters. Observationally, there are small tempera- 
ture inhomogeneities in the non-cooling flow clusters, ex- 
cept for merging clusters. This indicates that the heating 
source does not exist in all clusters, but exists only in the 
cooling-flow clusters. That is, the heating source may be 
closely related with the cooling properties. If clusters are 
self-similar within the cooling radius [equation (6)], the 
heating which accompanies the formation of cD galaxy's 
envelope is promising, since the existence of cD galaxies 
is strongly correlated with the existence of cooling flows. 

4. Discussion 

We considered the case that heating is an explanation of 
the "lack of cold gas" problem as an example. Of course, 
there are other possibilities. If inhomogeneous metallic- 
ity is the explanation, the cut-off temperature may be 
determined by the degree of inhomogeneity in the metal- 
licity, which may depend on the history of each cluster. 
Then, it is expected that the strong -/jcooi~?gas relation 
does not show up and the ratio Tcut/Tgas changes from 
cluster to cluster. We consider that observations from 
XMM-Newton can study these features. Next, we consider 
a case in which the absorption by cold (neutral) material 
is the explanation. In this case, the cut-off temperature is 
determined by the amount of cold gas, which may be pro- 
portional to the accretion rate, that is Lcooi- For brighter 
Lcooi clusters, the cut-off temperature may become higher. 
We think this tendency can be checked relatively easily in 
near-future observations. It is a future work to predict 
the relation between Lcooi and Tcut for other cases. 

In order to reproduce the observed 
luminosity-temperature relation, wc made the gas 
fraction vary with the cluster mass in the above dis- 
cussion. There is another possibility that the he gas 
profile deviates from the self-similar relation [equation 
(1)]. The self-similar relation [equation (1)] predicts 

1/2 

that Tc oc i?ci oc Tgas , which is supported by Vikhlinin 
et al. (1999) observationally; but this point is still in 
controversy. For example, Xu et al. (2001) showed 
that Tc oc Tgas and Mohr and Evrard (1997) insisted 
that the X-ray isophotal size is proportional to Tgas- 
These results indicate that non-gravitational processes 
are important in cluster evolution (e.g., pre- heating 
model, |00 [cite] cite. eh91 Evrard, Henry (1991); |00 
[cite]cite.kai91Kaiser (1991)). Deviations from the 
self-similar scaling relation [equation (1)] give useful 
information to understand cluster physics. If we consider 
these observed relations seriously, our discussion must be 
changed. Therefore, we performed the same procedure as 
that in sections 2 and 3, assuming Tc oc Rd oc Tgas instead 
of hq = constant. We obtained roughly the following 
relations: 

jio«A47^/Vrg-J, Tga.(xMy^ LboicxT|/2. (9) 
Thus, in order to reproduce the observed 



luminosity-temperature relation wc must abandon 
the assumption that the gas fraction is constant. The 
gas fraction must be proportional to Tgls oc M^Y*, then, 

2 74. 27 /8 

no oc Tgas ■ We thus obtain Lcooi oc Tgas , or steeper, 
depending on the first term in the brackets of equation (4) 
being much greater than 1 or not, instead of Tcooi oc Tj^^^ 
in section 3. If clusters are self-similar within the cooling 
radius [equation (6)], the heating rate, F, has to be 

proportional to M^^^^^, or steeper, instead of M^j in 
section 3. This difference may not significantly affect the 
rough discussion about the heating source in this Letter. 

On the other hand, recently, Ota (2001) suggested that 
clusters have two different length scales of about 50 kpc 
and 200 kpc for each core radius, assuming the Hubble 
constant to be 50 km s^^Mpc"^. Although we have not 
understood the origin of these two length scales, here, we 
tentatively consider that the larger core size describes the 
cluster properties. Since the core size does not correlate 
with the temperature strongly, we assume Tc to be con- 
stant as another example (it is hard to consider cluster 
size is constant, thus we do not assume Tc oc Rd) instead 
of no ~ constant. Then, we obtain roughly 

nocxA/^Y'aTgas, Rd'^Mll\ Tg,,^M^/\ (10) 

In this case, we obtain Tboi oc Tgis , thus the gas frac- 
tion must be proportional to Tgls oc M^-/^ in order to 
reproduce the observed luminosity-temperature relation. 

Then, Tip OC Tg/s and wc obtain Tcooi oc Tgas or steeper 
depending on the first term in the brackets of equation 
(4) being much greater than 1 or not. If clusters are self- 
similar within the cooling radius [equation (6)], the heat- 

7 /4 

ing rate has to be proportional to M^l , or steeper. This 
difference may not significantly affect the rough discussion 
about the heating source, too. 

Recently, Chandra observations have shown that some 
clusters have a complex X-ray surface brightness distri- 
bution, arising from the presence of radio lobes (e.g., the 
Perseus cluster: ^0 [cite] cite. fseOOFabian et al. (2000); 
the Hydra A cluster: OOO [cite] cite. mwnOOMcNamara et 
al. (2000); A2052: |gOO [cite] cite. bsmOlBlanton et al. 
(2001)). The size of the radio lobes is about several-times 
10 kpc, which corresponds to the size of the region where 
the "lack of cold gas" problem occurs. At present, there 
is no evidence for shock-heated gas surrounding the radio 
lobes. Therefore, the radio activity is not likely to be a 
major source of energy input. However, it is possible to 
affect our discussion through modifying the gas distribu- 
tion. We think that the connection between the effect of 
radio lobes on the surrounding gas and the "lack of cold 
gas" problem must be thought about in future work. 

At present, the "lack of cold gas" problem has been 
reported in only four clusters: A1795 (Tvir ~ 6 keV, 
000 [cite]cite.tam01Tamura et al. (2001)), Sersic 159-03 
(Tvir ~ 3 keV, 000 [cite] cite. kaaOlKaastra et al. (2001)), 
A1835 (Tvir ^ 8 keV, |00 [cite]cite.pet01Peterson et 
al. (2001)) and the Virgo cluster (Tvir ~ 2.5 keV, 000 
[cite]cite.mbflD01Matsushita ct al. (2001)). The contra- 
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dictory cases, i.e., cooling flow clusters without the cut- 
off temperature, are not yet reported by XMM-Newton. 
Thus, at present it may be premature to discuss the self- 
similarity; for example, studying the relation between Tcut 
and Tgas- However, it is worth pointing out that this 
problem has been reported in a relatively low tempera- 
ture cluster; Sersic 159-03, and the Virgo cluster. Thus, 
this feature must be universal and this point needs to be 
taken into account in considering the problem. If heating 
is the explanation, it is likely that the exponent of the 
heating, p, is around 2. 
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